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Abstract: By assuming that the scalar meson Aq (1430) belongs to the first excited states or the lowest lying ground 
states, we study the pure annihilation-type decays B —l A’q^( 1430)A^*^^ in the QCD factorization approach. Within 
the standard model, the branching fractions are at the order of 10“® —10“^, which is possible to be measured in the 
ongoing LHCb experiment or forthcoming Belle-II experiment. We also study these decays in the family non-universal 
Z' model. The results show that if mz' ~ OOOGeV — 0.02), both the branching fractions and CP asymmetries of 
B A'q"''( 1430)A'“ could be changed remarkably, which provides us a place for probing the effect of new physics. 


These results could be used to constrain the parameters 
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1 Introduction 

Heavy flavor physics has been a hot topic for many 
years in particle physics, since it is important for the 
study of CP violation and also a good place for searching 
of new physics signals. With more and more data from 
the LHCb experiment, many hadronic B decay modes 
are well studied experimentally, but not the case for the 
theoretical side. Among them, rare decays with flavor 
changing neutral currents are the most interesting, since 
they have a quite small branching ratios in the standard 
model, thus more sensitive to any new physics contribu¬ 
tions. 

In this work, we shall study the impact of a 
family non-universal leptophobic Z' boson on B ^ 
decays dominated by flavor changing 
neutral currents (FCNC). In the standard model, the 
FCNC processes are suppressed since it is only induced 
by loop diagrams. The branching fractions of these de¬ 
cays are predicted at the order of 10“® — 10“’^ [1], which 
is possible to be measured in the ongoing LHCb exper¬ 
iment or forthcoming Bellc-H experiment. While the 
family non-universal Z' boson leads to the tree-level fla¬ 
vor changing neutral currents, which may enhance the 
branching ratios of these decays. Motivated by this, 
many FCNC processes induced by Z' in flavor physics 
have been explored extensively in literatures [2-4]. Fur- 


of Z' model. 


thermore, as the scalar meson iFo^(1430) is now unclear 
whether it belongs to the first excited states (Scenario 1) 
or the ground states (Scenario 2) [5], we have to discuss 
both cases in the current work. 

Unlike the leptonic decays, the non-leptonic B decays 
are complicated because many degrees of freedom and 
scales are involved. Based on the effective field theories, 
there are three major QCD-inspired approaches for deal¬ 
ing with them, namely, the QCD factorization (QCDF) 
[6], perturbative QCD (PQCD) [7], and soft collinear 
effective theory [8]. In this work, we shall employ the 
QCDF for evaluating the relevant hadronic matrix ele¬ 
ments, as it is a systematic framework to calculate these 
matrix elements from QCD theory, which holds in the 
heavy quark limit ni}, —>■ oo. 

2 The family non-universal Z’ model 

The recent discovery at the LHC [9] of a light Higgs 
with mass around 125 GeV opened a new window to 
physics beyond the Standard Model (SM). In some new 
physics (NP) models [10], when the initial group breaks 
down to SU{2)i^ x U{1)y of SM, an extra group t/(l)' 
will usually be produced, which leads to an additional 
massive neutral gauge boson called Z'. If such a gauge 
boson were observed, it would be a concrete evidence 
of the existence of NP. Based on the assumption that 
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the Z' shares similar characters with Z boson of SM, 
many efforts have been made to search for Z' directly 
by analyzing the data of dilepton. Now, at the LHC, 
the lower mass limit is 2.86 TeV (1.90 TeV) at the 95% 
confidence level at 8 TeV colliding energy with an inte¬ 
grated luminosity of 19.5fb~^ by using e+e“ and 
[11] (or r+T“ [12]) events. However, such constraints 
from LHC are invalid if the Z' boson does not couple or 
couples very weakly with the leptons, thus one has to 
resort to hadronic channels. Although the couplings be¬ 
tween quarks and Z' are family universal in most models, 
the family non-universal Z' can also be realized in some 
models. For instance, such a family non-universal lepto- 
phobic Z' boson can be realized in Eq model [13]. The 
phenomenological studies of family non-universal Z' on 
possible colliders have been explored in detail recently 
[14]. 

On the gauge interaction basis, the interactions of the 
leptophobic Z' boson with SM quarks can be expressed 
as 


the effective Hamiltonian of SM [15], the operators of 
the forms {hd)v-A{qq)v-A and {bd)v-A{qq)v+A in eq.(3) 
have already existed in SM, so the Z' effect can be rep¬ 
resented by modifying the Wilson coefficients of the cor¬ 
responding operators. Thus, the eq.(3) can be rewritten 
as 

+h.c., (4) 

where Clan's 7 9 are the four quark operators existing in 
SM[15]. The additional contributions to the SM Wilson 
coefficients at the scale in terms of Z' parameters 
are thus given as 

=-3lfe , (5) 

AC9(r) = -(“) ' Bt, . ( 6 ) 


where the field ipi is the ith family fermion and Pl,r = 
( 1 =F 75)72 is the chirality projection operators, (e^ii) 
stand for the left-handed (right-handed) chiral couplings, 
and they are required to be hermitian because the La- 
grangian is real. When the weak eigenstates being ro¬ 
tated to the physical basis, the mass eigenstates will 
be obtained by 'ipL,R = fl'T'i r- Correspondingly, the 
coupling matrices of down-type quarks read 

B^ = V,, C, Vl , B2 = . ( 2 ) 

Since we do not need the couplings for up-type quarks, 
we will not discuss them here. Obviously, if the matrixes 
Cd,^ H are not proportional to the identity matrix, the 
nonzero off-diagonal elements in the B^’f appear, which 
will induce the FCNC interactions at the tree level. For 
simplicity, the right-handed couplings are often supposed 
to be flavor-diagonal. Then, the effective Hamiltonian of 
the b —>■ dqq transitions mediated by the Z' is given by 

V 2 Kgimz'J 

+By^{qq)v+A)+dc. , (3) 

where gi = e/(sindvv cosdvv) and niz' denotes the mass 
of Z' boson. The diagonal elements B^^^ are real due 
to the hermiticity of the effective Hamiltonian. In con¬ 
trast, the off-diagonal element B^^ might be a complex 
number with a new weak phase and such a newly 
introduced phase can be used in explaining the large di¬ 
rect CP asymmetries in B ^ Ktt [2, 3]. Compared with 


From above equations, we find that both the electro-weak 
penguins ACQ(r) and the QCD penguins ACs^s) will be 
affected by the new gauge boson Z'. Since the scale of 
new physics is expected to be much higher than that of 
electroweak scale, in order to show that the new physics 
is primarily manifest in the electroweak penguins, we 
therefore follow ref. [2-4] to assume ~ —2BPP. In 
this case, the Z' contributions to the Wilson coefficients 
are 

AC3(5) = 0 , AC9(r) = ^ ( 7 ) 

Ft6 Vtd 

where 


/ 32TOZ 


tdL jdX 

^bd^ss 

\gimz', 

I 

Vt,v:d 


dbd = Aig[Byy . 


{X = L,R), 


( 8 ) 

(9) 


Finally, we obtain the resulting effective Hamiltonian at 
the Mw scale 


njZ' 

^efF 


"Tr 

-l-h.c. 


/ g2mz 
\gimz' 


Bt 


■n 





Because all heavy degrees of freedom (including the Z') 
above the scale of W boson mass have been integrated 
out already and there is no new particles below niw, the 
renormalization group evolution of above new Wilson co¬ 
efficients down to low energies is exactly the same as in 
the SM [15]. 

Now, we will discuss the the ranges of new parame¬ 
ters and (j)bd- Because both the gauge group U{1)y 
and the new C/(I)' are expected to origininate from the 
same large group, the relation 32/111 ^ 1 is assumed. In 
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the experimental side, if the leptophobic Z' boson is de¬ 
tected at LHC, the mass of Z' should be about a few 
TeV, which means mzimz' is at the order of 0(10“^). 
In addition, the other parameters |i?^| and new 

weak phases (j)hd could be constrained by the data in¬ 
duced by FCNC. For instance, the mass difference B° 
and B^ requires ^ \VtbVf*J. Then, with experimen¬ 
tal data of Bd,s nonleptonic charmless decays, B^^^ ^ 1 
could be extracted. As for the new introduced phase 4>bd, 
recent discussions indicate that (pbd = —50° [3]. How to 
constrain these parameters globally is beyond the scope 
of this work. We will not discuss this topic explicitly 
here. So, to probe the new physics effect for maxi¬ 
mum range, we assume C = ^ [0.001,0.02], and 

(f’bd G [—180°, 180°], and set ( = 0.01 and cj)bd = —50° for 
the center values. 


3 QCD factorization calculation 

In the B —^ decays, none of the 

quarks in final states is the same as those of the initial B 
meson. Therefore they can occur only via annihilation 
type diagrams. However, in QCDF, the end-point singu¬ 
larity usually appears when calculating the annihilation 
type diagrams [6] . As a most popular way, the end-point 
divergent integral is treated as infrared sensitive contri¬ 
bution that is parameterized by 

where the unknown parameters Pa and (j)A could be fixed 
by the experimental data. This singularity can also be 
smeared by introducing the effective mass to gluon [16]. 
We will adopt eq.(ll) in this work. 

Within the framework of QCDF, the decay ampli¬ 
tudes of 5—^ {1A3Q)K^*'>^ can be written as 


A{B° ^K*+K-) 


A{B° K+K*-) 


A{B ^k;+k*-) = 


^ I /s/k* Ik [ {biSl + b4 + 64,ew)k*+k- + f ^4 - ^^4,ew^ ^ ] I , (12) 

^ /s/K*/if [(5 ii5]] + 54-I-54.ew)x+x*-+ f^4—2^4 ,ew^ ] k (13) 

^ -/s/kj/k* [(5iA[]-|-54-|-64^Bw)j^*+K*- + f^4-'r54,EW^ ji, (14) 

,= u.c I A ^ /k—K*+ J 


Gp 


A{B° -^K*+K^-) 


V2 




— fsfK^fK* [ (5lA[]+54 + ^4,Ew)j^.+ ^*- 


-h 



, ( 15 ) 


where the building blocks bi and bi^Ew read 

E, 

N2 

E, 

N? 


bi = 
1)2 = 
53 = 

bd = 
b3,EW = 
b4,EW = 


rCiAl, 

rC^Al, 


[CaAj-|-C5(A3-I-A3) + A^t-CgAg] , 

^[C'Al + C'Af], 

E [c'Hi+ c;(a;+A3O+jv^c'a;] , 

E[c[gHl + C'A']. (16) 


The expressions of the functions A];^ can be found in 
Ref.[5, 17], and the subscripts 1,2,3 denote the an- 


njhilation amplitudes induced from (V — A)(F — A), 

(F —A)(F-|-A) and {S — P){S + P) operators, while the 
superscripts i and / refer to gluon emission from the ini¬ 
tial and final-state quarks, respectively. When listing the 
two final state mesons M 4 M 2 in the formulas, we refer the 
first meson Mi to containing an anti-quark from the weak 
vertex and refer M 2 to containing a quark from the weak 
vertex. As for the aforementioned endpoint singularity 
Xa, we adopt the eq.(ll) with A = 0.5GeV. Note that 
the recent global fit of pa and to H —>■ SP,SV implies 
Pa = 0.15 and (f)A = 82° with = 8.3 [18]. Therefore, we 
set Pa G [0.1,0.2] and 4 >a £ [60°, 120°] in estimating the 
uncertainties. 

Because both B° and B could decay to 
Aro^(1430)Ar“ and A'o~(l'130)Ar+, we then define four 
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decay amplitudes, Af, Af, Af and Af as 

= {KS+K-m,Af = {KS-K+m; 

A, = (k;+k-, Af = {k;-k+ \b°). 


Then, the direct CP asymmetry is defined as 


\Af\^ + \Af\^-\Af\^-\Af\^ 

( 17 ) \Af\^ + \Af\^ + \Af\^ + \Af\^- 


(18) 


Table 1. Branching fractions of 5 —>• (1430)1^^*^ under different scenarios. 


Scenario 

Decay Modes 

SM 

Z' 


P"^Po*+(1430)iF- 

n q'7+0.43+0.23 
U.y 1 -0.31-0.12 

-1 Q7+0.57+0.35+0.62 
' -0.43-0.19-1.05 

Scenario-1 

W" ^k;-{uso)k+ 

Q qq~I”3.31-)-2.65 

0.00 2.55 i.39 

Q 1 7+3.25+2.60+0.59 
O-'L ‘ -2.50-1.36-0.18 

^k;+{U30)k*- 

n QQ+0.14+0.08 
U.OO o.10-0.04 

^ /ll^+0.16+0.10+0.10 

U.4:U o. 12-0.05-0.20 


W ^k;-{uso)k*+ 

14 ^74+5.75+4.54 
44_2.38 

1/1 /I 7+5.72+4.52+2.67 
i^.4: i -4.42-2.37-0.08 


b'" ^ k;+{U30)k- 

rj rO + 0.45+0.04 
U.OO o.29-0.02 

n fJC+0.50+0.08+0.07 
U.UO o.32-0.03-0.18 

Scenario-2 

^k;-{U30)k+ 

1 n7+0'^2+0.79 
-0.47-0.29 

-I pj^+0.71+0.77+0.08 
l.UO o.47-0.28-0.06 

W ^K;+iU30)K*- 

rj 1 1 +0.09+0.03 

U.ll o.06-0.02 

rj 1 9+0.09+0.04+0.01 

U.iZ-0.06-0.02-0.02 


W ^k;-{uso)k*+ 

1 0/1+1-26+1.38 
1.0^0.83-0.52 

n.+1.25+1.37+0.04 
l.O^0.82-0.51-0.04 


In the experimental side, the four time-dependent de¬ 
cay widths are given by {f = K~ and f = Kq~K^) 




[1 -I- cos Amt — Sf sin Amt], 

[1 — Cjcos Amt -I- S'/ sin Amt], 

e-"‘i(|A/f+ |i/f) 

[1 -I- CjCOS Amt — S/ sin Amt], 

e-"‘i(lA/r + li/n 

[1 — C/ cos Amt + Sf sin Amt] .(19) 


In above functions. Am means the mass difference of 
B°/B meson, and T is the average decay width of the 
B meson. The auxiliary parameters Cf and S/, that can 
be extracted from the data, are given by 


Gf = 

l+l^ + l+P’ 

(20) 


2Im(A/) 

(21) 

— 

1 + 4//4/ 2’ 

= 

Af 

(22) 

V:,Vta Af- 


The definitions of C/ and S/ can also been obtained by 
replacing / with /. In order to show the implications of 
above four parameters, we usually use the following four 


new parameters: 


G — f^{Cf + Cf), 

AC = \{Cf-Cf), 

(23) 

S = \{Sf + Sf), 

AS=\{Sf-Sf). 

(24) 


Physically, S and C are related to the mixing-induced 
CP asymmetry and the direct CP asymmetry, respec¬ 
tively. Moreover, AC and AS are CP-even under CP 
transformation A/ —>■ 1/A/. 

4 Numerical Results 

Using the parameters of refs. [5, 17, 18], with and 
without Z' boson, we present the predicted branching 
ratios of P —?► iCo*(1430)iC^*^^ under different scenar¬ 
ios (Pro^(1430) as the first excited states or the ground 
states) in Table. 1. The predictions of SM are also listed 
for comparison. Within the PQCD approach based on 
the kr factorization, Liu Xin et.al. had studied these 
decay modes [19]. Comparing their results with ours, 
we find that their results of branching ratios are much 
larger than ours by 1 ~ 2 orders of magnitude. Currently, 
we cannot determine which approach is better, but we 
hope the future experiment can test these two different 
approaches. 

From Table. 1, for decay modes B —>■ 

/Co + (1430)iF^*U^ the differences between SI and S2 is 
very small, so it is hard for us to discriminate two dif¬ 
ferent scenarios using these two decays. In contrast, the 
branching fractions of —>■ Kq~{ 14:30)KiS+ under SI 
and S2 have large differences, so if these two modes can 
be measured precisely they may be used to determine 
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whether A"o(1430) belongs to the ground states or the 
first-excited states. 

In Table. 2, we also list our predictions of Acp, C, 
AC, S and AS for the final states Kq^K~ and Kq^K*~, 
under two different scenarios in both SM (the upper val¬ 
ues) and Z' model (the lower values). For all results, the 
first errors are from the uncertainties of decay constants 
and the light-cone distribution amplitudes of final states 
and the second errors come from the Pa and (1 )^, and the 
last errors in the Z' model are the results by scanning 


the possible ranges of C and (j)bd- From the numerical re¬ 
sults, we find that the largest uncertainties are from the 
Pa and /pA- This is in contrast with other decay modes 
dominated by the spectator diagrams, whose uncertain¬ 
ties taken by the two above parameters are small. Unlike 
branching fractions, the CP asymmetry parameters are 
not sensitive to non-perturbative hadronic parameters, 
where these uncertainties are canceled because they are 
ratios of different amplitudes. Therefore, these parame¬ 
ters are more sensitive to the effect of NP. 


Table 2. The CP-violating parameters Acp, C, AC, S and AS (%) of B ^ Kq{U30)K^*\ 


Decay 

s. 

Acp 

c 

AC 

s 

AS 

-H 

SI 

r\ 7Q+0.01+0.01 

u. < y_o.oi-o.oi 

rj 71 +0.01+0.01+0.22 

U. ( -L-0.01-0.01-0.11 

A A/I+0.01+0.02 

U.U+_0.01-0.01 

A AA+0.00+0.00+0.05 

U.UZ o.00-0.01-0.11 

A AA+0.01+0.02 
U.UU_o.02-0.03 

A Ai +0.01+0.02+0.22 

U.U 1-0.01-0.03-0.12 

A A/I+0.02+0.01 

U.U+_o.02-0.01 

A A/I+0.02+0.01+0.02 

U.U+o.02-0.01-0.07 

A n/i+0.03+0.02 

U.0+_o.03-0.02 

A QO+0.02+0.01+0.19 
U.OO o.02-0.02-0.03 

S2 

A nQ+0.06+0.20 
U.Zy_o.o9_0.14 

A Q/1+0.07+0.18+0.14 
U.Z+o.07-0.12-0.05 

A aa+0.04+0.04 
U.UZ_o.o2-0.05 

A AA+0.03+0.04+0.08 
U.UZ o.01-0.05-0.19 

n oc:+0.04+0.09 
U.ZJ_o.o2-0.12 

A 97+0.06+0.12+0.10 

U.Z/ -0.05-0.15-0.16 

A q/::;+0.05+0.06 
U.OU_o.05-0.04 

A q c+0.06+0.06+0.08 

U.OO o.06-0.04-0.21 

A 9A+O.I2+O.IO 

U.ZU_o.iO-0.16 

A 9A+0.14+0.11+0.19 

U.ZU o. 11-0.18-0.05 

a:*±a:*=f 

SI 

A AA+0.00+0.00 

U.i^U_0.00-0.00 

A 0/1+0.00+0.00+0.03 

U.y+-0.00-0.00-0.01 

A Ai +0.00+0.01 

U.Ul-o.oO-0.00 

A AI +0.00+0.01+0.04 

U.U1 0.00-0.00-0.02 

A A9 + O.OO+O.OI 

U.UZ_0.00-0.01 

A AI +0.00+0.01+0.07 
U.U 1-0.00-0.01-0.04 

A Ai +0.00+0.00 

U.Ul-o.oO-0.00 

A AI +0.00+0.00+0.01 

U.Ul-0.00-0.00-0.00 

A aq+O.OO+0.02 

U.UO_0.00-0.00 

A QO+0.00+0.00+0.02 

u.yo o.00-0.00-0.01 

S2 

A qo+0.02+0.03 
U.OO_o.oi_0.02 

A Q7+0.02+0.03+0.02 

U.O/ -0.01-0.02-0.01 

A A1 +0.01+0.02 

U.Ul_o.oi_0.02 

A AQ+0.02+0.03+0.02 

U.UO o.01-0.02-0.06 

A 1 1 +0.03+0.06 
U.11_o.o2-0.06 

A 1/1+0.04+0.07+0.02 
U.I+-0.03-0.07-0.08 

n A9+O.OI+O.O8 
U.UZ_0.01-0.04 

A A9+0.01+0.01+0.02 

U.UZ o.02-0.01-0.02 

A 01 +0.02+0.02 

U.i^l_0.02-0.02 

A 09+0.03+0.03+0.02 

U.yz 0.02-0.03-0.02 


In order to show the effects of Z' clearly, we also plot 
the variations of the branching ratios as functions of the 
new weak phase pbd with different (p = 0.001,0.01,0.02 
under different scenarios in Fig.l. Including the newly 
introduced Z' boson, one can see that if C < 0-01 the 
effects of Z' are not large enough to be detected in these 
four decay modes, because the new physics contribu¬ 
tions are buried by the large uncertainties of hadronic 
parameters. If C is around 0.2, the branching fraction of 

—>■ Arg^(1430)iF“ under SI will be changed remark¬ 
ably to reach to 2.0 x 10“’', which could be measured in 
the forthcoming Belle-II experiment. 

The relations of the direct CP asymmetries A^p, C, 
AC, S and AS” with with different ^ are also pre¬ 
sented in Fig. 2 and Fig. 3, for B —>■ A'g*(1430)A"^ and 
B —7> Aro"‘"(1430)Ar*"F, respectively. From Fig. 2, the ob¬ 
servables oi CP asymmetries are much sensitive to the 
Z' than the branching fractions. For example, under SI, 
the direct asymmetry of ^ ArQ*(1430)Ar^ is 79% in 
SM, while it would change to 49% with Z' boson. In fu¬ 
ture, these observables could be used to probe the effect 


of new physics. If the Z' were detected in the colliders 
directly, these decays would also be useful to constrain 
the couplings. 

5 Summary 

Within the QCD factorization approach, we have 
studied the pure annihilation type decays B —>■ 

A:*±(1430)A:(d=F 

in SM and the family non-universal lep- 
tophobic Z' model. Both the branching fractions and the 
CP asymmetry observables have been calculated. The 
branching fractions we predicted are 1-2 orders of mag¬ 
nitude smaller than the results from PQCD approach. 
When the Z' involved, if ruz' > ITeV (C < 0.01), its 
contributions will be buried by the large uncertainties of 
SM. If iTiz' ~ OOOGeV (C = 0.02), both the branching 
fractions and CP asymmetries of B° —>• Aro'''(1430)Ar“ 
could be changed remarkably, which provides us a place 
for probing the effect of new physics. These results are 
hopeful to be tested in Belle-II, LHC-b or the future high 
energy circular colliders. 
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Fig. 1. Under different scenarios(Sl and S2), the branching ratios of B ^ A'q^( 1430)A'^*^^ as functions of the 
weak phase (fibd, the dotted (red), dashed (orange) and dotdashed (blue) lines represent results from the ^ = 
0.001, 0.01, 0.02, and the solid lines (black) are the predictions of SM. 
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Fig. 2. Under different scenarios(Sl and S2), the CP-violating parameters Acp, C, AC, S and AS {%) of B ^ 
isro(1430)A' as functions of the weak phase (f)bd, the dotted (red), dashed (orange) and dot-dashed (blue) lines 
represent results from the ^ = 0.001, 0.01, 0.02, and the solid lines (black) are the predictions of SM. 
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Fig. 3. Under different scenarios (SI and S2), the CP-violating parameters Acp, C, AC, S and AS (%) of B ^ 
Kq{1A:?>0)K* as functions of the weak phase <j)bd, the dotted (red), dashed (orange) and dot-dashed (blue) lines 
represent results from the C = 0.001, 0.01, 0.02, and the solid lines (black) are the predictions of SM. 
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